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Abstract 

We construct a TV— dimensional Gaussian landscape with multiscale, trans- 
lation invariant, logarithmic correlations and investigate the statistical me- 
chanics of a single particle in this environment. In the limit of high dimension 
N — > oo the free energy of the system and overlap function are calculated ex- 
actly using the replica trick and Parisi's hierarchical ansatz. In the thermody- 
namic limit, we recover the most general version of the Derrida's Generalized 
Random Energy Model (GREM). The low-temperature behaviour depends 
essentially on the spectrum of length scales involved in the construction of 
the landscape. If the latter consists of K discrete values, the system is char- 
acterized by a K-step Replica Symmetry Breaking solution. We argue that 
our construction is in fact valid in any finite spatial dimensions N > 1. We 
discuss the implications of our results for the singularity spectrum describ- 
ing multifractality of the associated Boltzmann-Gibbs measure. Finally we 
discuss several generalisations and open problems, such as the dynamics in 
such a landscape and the construction of a Generalized Multifractal Random 
Walk. 

PACS numbers 64.60.De, 64.60.A1 
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1 Introduction 



Ever since the seminal paper of Goldstein in 1969 [1], the idea of energy landscapes per- 
vades the theoretical description of glasses, disordered systems, proteins, etc., see [2, 3, 4, 5] 
and references therein. The general idea is to describe the statics and dynamics of the 
whole system, or one of its subparts, by a single point particle moving in a random po- 
tential, which encodes the complexity of the original system. The hope then is to be 
able to classify the possible classes of random potential and to establish generic, universal 
properties, in the spirit of Random Matrix Theory. In this respect, the Parisi solution for 
spin-glasses is fascinating: it reveals that in this case the energy landscape has a surpris- 
ingly complex, hierarchical structure of valleys within valleys within valleys, etc. [6]. It is 
often argued that this construction is very specific not only to infinite range spin models, 
but also to infinite dimensional landscape models. In particular, the ultrametric proper- 
ties of Parisi landscapes seem at first sight hardly compatible with a finite dimensional, 
translation invariant random function. 1 In this paper, we provide an explicit construction 
of a Gaussian random potential in Euclidean, N dimensional spaces, with a specific form 
of long-ranged correlations which reproduces all the features of Parisi landscapes. More 
precisely, we show that the thermodynamics of a single particle in a multiscale, logarith- 
mically correlated potential is exactly described by Derrida's Generalized Random Energy 
Model (GREM, [8]), with an arbitrary (possibly infinite) number of levels of hierarchy. Al- 
though our proof concerns, strictly speaking, the limit N — ► oo, we are confident that our 
results hold in arbitrary finite dimension N > 1. This conviction is built both on physical 
arguments and on the beautiful results of Carpentier and Le Doussal [9] on the monoscale 
version of our model in finite dimensions, which, as shown recently, match precisely the 
exact results of the same model when N — > oo [10]. 

As is well-known, the Gibbs-Boltzmann measures in systems with disorder often pos- 
sess the interesting property of being multifractal, see the papers [9, 11, 12]. In fact, 
this property is not unrelated to the multifractality of the wavefunctions in disordered 
electronic systems (see [13] for a comprehensive discussion of the last topic and further 
references). The investigation of multifractal measures of diverse origin has been a very 
active field of research in various branches of physics for about two decades now[15, 14, 13]. 
From this point, we show that our results imply, in particular, the possibility of a rather 
rich and unusual behaviour of the singularity spectrum describing multifractality of the 
Gibbs-Boltzmann measure arising in our multiscale logarithmic model. 

Another closely related aspect is that the monoscale model of Carpentier and Le Dous- 
sal is known to be the building block in the construction of an exact multifractal random 
walk (MRW), proposed by Bacry, Muzy & Delour to describe financial time series [14]; the 
Boltzmann weight in one language corresponds to the local volatility in the other. So our 

^dn fact, random potentials with a hierarchical Parisi structure can be constructed in finite 
dimensional space in a kind of ad-hoc way, by following step by step the real-space interpretation 
Replica Symmetry Breaking: see [7] for a discussion of this point. 



2 



extended multiscale model can also be interpreted as the construction of a Generalized 
Multifractal Random Walk (GMRW), in the same sense as the GREM generalizes the 
Random Energy Model. Physically, the result is that the n-th moment of the distance 
travelled by the random walk scales with an exponent that does not only depend on n but 
also on the epoch, ie the logarithm of the time lag. 

The outline of the paper is as follows: in Sect. 2 we introduce landscape models in 
full generality, recall the set of previously established results and discuss how they can 
be understood qualitatively. We then focus on the case of logarithmically correlated 
landscapes, summarize the recent findings of ref. [10] and their relation with the results 
of Carpentier and Le Doussal [9]. In Sect. 3, we define precisely our multiscale random 
landscape model, discuss its physical motivation and detail our analytical calculations in 
the large dimensional N — > oo limit, where we recover exactly the GREM results. We end 
this section by describing the most interesting features of the multifractality spectrum 
of the associated Boltzmann-Gibbs measure implied by our results for the multiscaled 
landscapes. Finally, in Sect. 4 we put forward several conjectures about the dynamics 
of a point particle in such landscapes, open problems and generalisations, concerning 
for example Generalized Multifractal Random Walks. Some more technical points are 
relegated to the Appendix. 

A short account of some of the results of the present paper was presented in [16]. 

2 Thermodynamics of a particle in a random po- 
tential 

2.1 General discussion 

As alluded to in the introduction, the "toy-model" of a classical particle in a random 
potential exhibits a rich variety of behaviour which mimics many of the dynamical and 
thermodynamical properties of glassy systems. The model is defined as follows: the po- 
sition of the particle, confined inside an A— dimensional spherical box of radius L, is 
described by the coordinate vector r = (n, rjv), |r| < L. It feels a random potential 
V(r), which we conventionally choose to be Gaussian-distributed with zero mean, and 
with covariance chosen to be isotropic, translation invariant and with a well-defined large 



In Eq.(l) and henceforth the notation {■ ■ -) v stands for an ensemble average over the 
random potential, and f(u) is a function of order unity belonging to the so-called class 
described in detail, e.g., in the book by Yaglom [17]. The functions f(u) € 2?oo are 
such that they represent covariances of an isotropic random field for any spatial dimen- 
sion N > 1. There are two essentially different types of such functions. The first type 



A-limit: 
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corresponds to genuine isotropic random fields, and those f(u) are characterized by a non- 
negative normalizable "spectral density function" f(k) > 0, k > in terms of which f(u) 
is represented as (see [17] p. 354): 

f(u) = [°° e- k2u f(k)dk, /(0) = [°° f(k)dk < oo . (2) 

JO JO 

In particular, f(u) is decreasing and convex, i.e. satisfies f'(u) < 0,f"(u) > \/u > , 
and in addition f'(u — > oo) = 0. Here and below the number of dashes indicates the order 
of derivatives taken. A few important families of such functions listed in [17] are, e.g., (i) 
f(u) = Ce~ au \ < 7 < 1, C >0, a > (ii) /(it) = C/(a + u) 7 , 7 > 0, C > 0, a > and 
(hi) f(u) = C{au) 1 / 2 K 1 (y/au), 7 > 0, C > 0, a > 0, where if 7 (x) stands for the modified 
Bessel (a.k.a. Macdonald) function. In the physical literature the random fields of that 
type are frequently called potentials with short-ranged (SR) correlations. 

The second type of covariances occurs in the situation when the normalization integral 
Jo°° /(fe)dfc diverges. It corresponds to long-ranged (LR) random fields with isotropic 
increments also known as locally isotropic random fields (see e.g. [17], p. 438). The spectral 
density function now must satisfy the condition / °° k ^ +l f(k)dk < 00 which allows one to 
prove that in any dimension N > 1 there exists a random field whose structure function 
i((V(0) - V{v)f) v = /(0) - f(u) is given by 

f°° ,2 - 

f(0)-f(u)= dk{l-e~ k U )f{k)dk + Au, A>0. (3) 

JO 

In what follows we will impose an additional requirement f'(u — ► 00) = 0, which ensures 
A = (no external driving force acting on the particle) . It is also easy to convince oneself 
that in the present model the difference between the covariance and the structure function, 
i.e. the value of /(0), is immaterial for the free-energy calculations. The most widely- 
known example of the locally isotropic LR field is the so-called self-similar random field, 
see [17] p. 441, characterized by the spectral density f(k > 0) = fc -27-1 , < 7 < 1. The 
corresponding covariance behaves as 

f{u) = f{Q)-C 1 u<. (4) 

In particular, for N = 1 and 7=1/2 this is the example of a simple Brownian motion for 
the potential, corresponding the celebrated Sinai model. 

After specifying in detail the class of random potentials involved in our construction, let 
us turn to thermodynamics of the model characterized by the following partition function 
and the corresponding free-energy: 

F = -T(\nZ) v , Z=[ exp -f3V(r) dr , (5) 

J\r\<L 

where (3 = y stands for the inverse temperature. A variant of the model consists in 
replacing the spherical box |r| < L by a confining harmonic potential — /ir 2 /2. Such a 
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model has been studied extensively since the mid-eighties. It was originally proposed 
in N = 1 dimension as a toy-model for a randomly pinned domain wall [18] or of a 
directed polymer in a random potential [19] and studied using a variety of methods [20, 
21, 9, 22], some of them being exact. Another case where analytical calculations can be 
performed is the high-dimensional limit N — > oo [23, 24], where a Gaussian variational 
ansatz with Replica Symmetry Breaking becomes exact. One finds that the nature of the 
low temperature phase is essentially dependent on the behaviour of the covariance at large 
distances [23, 24]. Namely, for typical short-ranged correlated potential the description of 
the low temperature phase was found to require the so-called one-step replica symmetry 
breaking (1RSB) scheme of Parisi. In contrast, for the potentials growing as u 7 , see Eq.(4) 
, the full infinite-hierarchy replica symmetry breaking (FRSB) scheme has to be used. In 
fact, this problem was reconsidered recently in [10] in much detail using an alternative 
method that directly focusses on the degrees of freedom relevant in the limit N — > oo, and 
employs the Laplace (aka saddle-point) method for evaluating the integrals. In the limit 
N — > oo, one actually finds a true phase transition as a function of temperature provided 
the radius of the confining sphere L is scaled as Rs/N. The effective size R < oo (which is 
accidentally just half of the length of an edge of the cube inscribed in this sphere) is then 
used as the main control parameter of the model. The chosen scaling L ~ y/N formally 
stems from the property of the argument of the correlation function Eq.(l) to become of 
order of unity for separations of order of y/~N. More importantly, it simultaneously ensures 
that the volume Vl = tt n / 2 y(n/2+i) °^ our spherical sample retains in the limit N 3> 1 
the natural scaling with size R and dimension iV: In Vl = Nln (R/Rq)+ smaller terms, 
with Rq being a constant of order of unity. Such a behaviour is essential since the phase 
transition is physically induced by a competition between entropic effects, which tend to 
delocalize the particle over the sphere, and the minima of the random potential which 
tend to attract, and possibly to localize the particle over a finite number of favorable sites. 
In the short-range case, the number of effectively independent "sites" is of the order of 
Vl oc R n which ensures, thanks to the N factor in front of f(u) in Eq. (1), that the 
minimum of the Gaussian potential scales as s/N Vln R N oc N. This indeed can compete 
with the entropy of the order N In R. All these arguments demonstrate that indeed R is 
the most natural measure of the sample size. 

The fact that a true thermodynamic transition exists for a finite sample size R < oo 
is obviously a somewhat pathological feature of the limit of infinite dimension N — ► oo 
taken first. Indeed, in that limit the total number of thermodynamic degrees of freedom is 
infinite even for finite R. At the same time, at any finite spatial dimension N < oo phase 
transitions only may occur in the thermodynamic limit of infinite sample size R — > oo. 
From this point of view it is natural to inspect the R — ► oo behaviour of the transition 
temperature T C (R) at which the system experiences a continuous breakdown of replica 
symmetric solution. For the present model it reads, according to [10]: 




(6) 
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in agreement with a similar result for the confining quadratic potential case [23, 24]. 
Therefore, one finds that T C (R) tends for R — ► oo either to zero in the SR case 2 , or 
to infinity in the LR case. A more detailed analysis shows correspondingly that the 
large R behaviour of the free-energy is FiT^R^^ ~ -T\nR N in the SR case where 
entropy dominates, and F(T)|^_> ~ — NRy in the LR case where the deepest minimum 
dominates. 

2.2 Logarithmically correlated potentials 

From Eq. (6) above, T c (i?)|^oo appears to have a well defined limit when f"(u) ~ u~ 2 
for large u, corresponding to a logarithmically growing correlation function of random 
potential, a case overlooked in previous studies [23, 24]. The peculiarities of that case 
can be traced in a few different ways. To this end it is appropriate to mention a precise 
mathematical criterion proposed recently in [10] to classify statistical mechanics behaviour 
induced by SR vs LR correlated potentials. The criterion uses the notion of the so-called 
Schwarzian derivative {f'(u),u} = —S(u)/[f"(u)] 2 , where S(u) is expressed in terms of 
f(u) as 

S(u) = l[f'"(u)] 2 -f"(u)f""(u). (7) 
In terms of S(u) it was demonstrated that 

• any potential whose covariance function satisfies the condition 

S(u) > V-u > must have a 1RSB low temperature phase. It is easy to check that 
such a situation includes, in particular, the standard families of the SR potentials (i) 
and (ii) listed after Eq.(2). For a general case of SR fields with covariances defined 
via Eq.(2) one always has S(u) > for large enough values of u, which is the most 
essential range in the thermodynamic limit. 

• Any potential whose correlation function satisfies the condition 

S(u) < Vu > must necessarily have the FRSB low temperature phase. This 
condition holds for the standard LR correlation functions of the type Eq.(4), i.e. 
for f(u) = f - g 2 (u + a 2 ) 7 , with < 7 < 1 and / - g 2 a 2 "' > 0. It is natural to 
conjecture that typical LR random fields with independent increments should be of 
this type for large enough u. 

Clearly, the above criterion naturally singles out as a special marginal case random 
potentials satisfying S(u) = 0. The only function satisfying this condition globally, i.e 
for all u > 0, and satisfying also the requirement f'(u — > 00) = is indeed given by a 

2 The situation is slightly more complicated, as for SR case there exists another temperature 
Ti(R) > T C (R) where the system experiences a discontinuous breakdown of replica symmetry. 
However one can show that Ti(R — > 00) — > as well [10]. 
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logarithmic correlation function, of the type considered by Carpentier and Le Doussal in 
finite dimensions [9]: 

f(u)=f -g 2 ln(u + a 2 ), (8) 

where g and a are given constants, and fo is such that /(0) = fo — g 2 In (a 2 ) > 3 . Let us 
stress that the expression Eq.(8) is a legitimate covariance function belonging to the 
class of LR locally isotropic fields, Eq.(3). Indeed, it corresponds to the spectral density 
of the form f(k) = ^f-e _a2fc2 , which satisfies the required condition f °° k ^ +l f(k)dk < oo. 

In such a case of logarithmically-correlated potentials the solution found in [10] has 
features of both the SR-1RSB and LR-FRSB regimes. The critical temperature T c remains 
finite for large systems, and is given by: 

T C (R ^oo)=g. (9) 

Physically, the minima of the potential now typically behave as: 



V min (R) ~ -gy/N\/2]nR N ^/ln(a 2 + R 2 ) ~ -2gN In R, (10) 

while the entropy contribution is —TNlnR, suggesting that indeed some change of physics 
should take place when T ~ g (sec [9] for further elaboration of this argument). 

Furthermore, it is easy to check that the free-energy expression found in [10] for arbi- 
trary R and a reduces in the limit R » a to that of the famous Random Energy Model 
[27]: 

1 Fm , / T(l + g 2 /T 2 )lnR, T > T c 

-^F(T)\R»a - j 2glnR ^ T <T C (U) 

Interestingly, these results coincides precisely with the Renormalisation Group results of 
Carpentier and Le Doussal in finite dimensions (up to a rescaling of their coupling constant 
g = sjo by a factor y/~N, as indicated by Eq. (1)). The interpretation is the same as for the 
REM: below T c , the partition function becomes dominated by a finite number of sites where 
the random potential is particularly low, and where the particle ends up spending most of 
its time [25, 26]. For a more quantitative description of the particle localization, useful in 
the following, it is natural to employ the overlap function defined as the mean probability 
for two independent particles placed in the same random potential to end up at a given 
distance to each other. Denoting the scaled Euclidean distance (squared) between the 
two points in the sample as V, and employing the Boltzmann-Gibbs equilibrium measure 
Pp(r) = zy3j ex P — f3V( r ) the above probability in thermodynamic equilibrium should be 
given by 

*<P) = ll dr lP/3 (n) / dr 2 ^(r 2 H(p-^| ri -r 2 | 2 )\ (12) 
W|ri|<L J\r 2 \<L \ // 



3 As noted above, the value of /(0) is irrelevant for the thermodynamics of the system, so that 
.fo can be dropped from the calculations. We systematically disregard such constants in the rest 
of the paper. 
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where here and henceforth 5 denotes the Dirac's 5-function. The disorder averaging in 
(12) can be calculated following the same standard steps of the replica approach as the 
free energy itself. For convenience of the reader we sketch the procedure for the present 
model in the Appendix A. 

With the function tt(V) in hand we can ask, in particular, what is the probability for the 
particle in logarithmically correlated potential to end up at V = 0(a 2 ), i.e. at a distance 
of order of the small cutoff scale. The answer turns out to be zero in the high-temperature 
phase T > T c , confirming the particle derealization over the sample. In contrast, in the 
lower temperature phase T < T c the probability is finite: it (0(a 2 )) = 1 — T/T c , again in 
full agreement with REM calculation [27, 28, 29, 9]. 

A logarithmic growth of the variance of the potential might look an academic oddity, 
but in fact is not, and appears naturally in various systems of actual physical interest. We 
warmly recommend the paper of Carpentier and Le Doussal [9], which discusses in detail 
the connection to many other interesting and important physical problems, like directed 
polymers in random environment [30], or a quantum particle in a random magnetic field 
[11]. In the next section, we introduce and study a very natural, multiscale generalisation 
of this model. 

3 A multiscale logarithmic potential 

3.1 Motivation and definition of the model 

The main observation of the present paper is that the above picture, despite looking rather 
complete, still misses a rich class of possible behaviour that survives in the thermodynamic 
limit R — > oo. Namely, given any increasing positive function &(y) for < y < 1, we 
demonstrate below that if one considers potential correlation functions f(u) which take 
the following scaling form 



the thermodynamics of our system in the limit R — ► oo is precisely equivalent to that 
of celebrated Derrida's Generalized Random Energy Model (GREM)[8]. The REM-like 
case discussed above turns out to be only a (rather marginal) representative of this class 
corresponding to specific choice of the scaling function $(?/) = g 2 y. 

Let us explain the motivation of the above form, which will make the physical inter- 
pretation of the results (as well as some technical calculations) quite transparent. The 
idea is to write V(r) as a (possibly infinite) sum of K independent Gaussian potentials: 




(13) 



K 




(14) 



i=i 
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each with a simple logarithmic covariance as in (8): 

(V i (r 1 )V j (r 2 )) v = 6 i , j Nf i ^(r 1 -r 2 ) 2 ^ , = -g 2 In (u + a 2 + a 2 ) , (15) 

each with its own strength constant gi and small-scale cutoff a^, which we choose to grow 
as a power-law of the system size: 4 ai = R Vi with < V{ < 1. Taking the continuum limit 
K — ► oo with a certain density p(z/) of exponents fj, we end up with: 

f(u) = - j\{v)g 2 {v)\n(u + a 2 + R 2v ^jdv, 0<x<R 2 . (16) 

Now, introducing u + a 2 = R 2y and identifying with Eq. (13) in the R — > oo limit, we find 
that the function $ has the following representation: 

Hy) = V I" p{v)g 2 {v)dv+ f up(u)g 2 (u)du, (17) 

JO Jy 

the previous REM case corresponding to p(y) = 5{y). Note also that in this representa- 
tion, <3?'(y) = Jq p{v)g 2 (v) dv > 0, and also <&"(y) > 0. The main result of this work is the 
following: depending on the nature of the spectrum of the exponents is, discrete or con- 
tinuous, we will recover, in the thermodynamic limit, either the free energy of the original 
GREM with discrete hierarchical structure, or of its continuous hierarchy analogue (see 
(32) below) analysed recently in much detail by Bovier and Kurkova [31], and appearing 
also in earlier studies of random heteropolymers[32]. 

The physical interpretation of our results also generalize the discussion of the previous 
section, Sect. 2.2, in a natural and, we believe, rather beautiful way. Instead of one 
localisation transition temperature T c where the particle chooses a finite number of "blobs" 
of size O(a) where the potential is particularly deep, there appears K different transition 
temperatures, where the particle localizes on finer and finer length-scales. The largest 
transition temperature T\ corresponds to a condensation of the Boltzmann-Gibbs weight 
inside a few blobs of large size O(R), but the particle is still completely delocalized inside 
each blob. As the temperature is reduced, the REM condensation takes place over smaller 
blobs of size 0(R U ) inside each already occupied large blobs, and this scenario repeats 
itself as the temperature is reduced, each time "zooming" in on a smaller scale 5 . To see 
this most clearly we quote the simplest example beyond REM, the two-scale logarithmic 
model characterized by the density of exponents g 2 (v)p(v) = g 2 5(v) + g 2 5(u — ui), with 
< v\ < 1. The system turns out to be described by two different critical temperatures 

T\ = \J g\+ g 2 > T 2 = g 2 - As we will be able to show, see Eq.(45), in our model the 

4 One could in fact multiply this power law behaviour by a slow function of R with no impact 
on the following results in the limit mi? — > oo. 

5 see [45] for a related discussion of the idea that temperature plays the role of a microscope in 
the context of spin-glasses. 
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probability for two particles to be found at a distance V = 0{R 2u ) apart is given by: 



[ 5{y - 1) , T > Ti 

n(i/) = | + T 2 <r<r 1 , (i 8 ) 

( (l-g)«5( l /) + (^- ^^-^O + ^-l), 0<T<T 2 

The first two lines reproduce the former REM scenario, with Ti standing for T c and the 
scale ai = 0(R U1 ) playing for T > T 2 the role of the lowest discernible cutoff scale. The 
last line describes quantitatively the "zooming in" from the scale a\ = 0(R Ul ) to the even 
smaller cutoff scale a = O(R ), which becomes discernible below T = T 2 and dominates 
more and more when T — > 0. 



3.2 Analytical results for N — > oo 



We aim to compute the equilibrium free energy per degree of freedom of our model, i 7 ^ = 
lhn/v— >oo Fn/N, where Tjy is defined in Eq. (5). The disorder average is performed in a 
standard way using the replica trick. The replicated partition function {Z n )y is evaluated 
exactly for \ < n < N — 1 m. the large- N limit by the Laplace method, after exploiting a 
high symmetry of the integrand stemming from the symmetry of the correlation function 
Eq.(l). The replica limit n — > is then performed in the standard framework of the Parisi 
hierarchical ansatz. The details of the corresponding analysis can be found in [10], and 
we give below a summary of the most essential formulae. 

From the point of view of the Schwarzian derivative criterion recalled above, our models 
(13) and (16) are such that for any finite R < oo, the low-temperature phase is charac- 
terized by continuous FRSB Parisi pattern with infinite level of hierarchy. This holds 
invariably, even when we make the choice of a discrete set of K > 2 distinct exponents 
< vk < ■ ■ ■ < v\ < 1. Indeed, the function S(u) defined in Eq.(7) when calculated from 
Eq.(16) reads: 



S(u) 



(u+slfiu + slY 



1 2 



{u + sl) {u + S„i) 



dv dv\ 



(19) 



where we have used the short-hand notation s. 



= a 



+ R 2u . The expression in the right- 



hand side of (19) is manifestly strictly negative, apart from the discrete K = 1 case when 
it is zero. Only in the thermodynamic limit R — > oo shall we find that S(u) — ► for most 
values of u. In this limit the system effectively recovers a GREM structure corresponding 
to the replica symmetry breaking pattern with K levels of hierarchy 6 . 

For finite R, the low temperature phase is therefore characterised by the existence 
of a nontrivial, non-decreasing function x(q), q G [qo,qk], with the two parameters go 
and qk satisfying the inequality < go < qu < Qd = R 2 - The corresponding can 



6 An alternative construction of multiscale logarithmic landscapes which show replica symmetry 
breaking with exactly K levels of Parisi hierarchy for finite R is proposed in the Appendix B. 
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be written in terms of only those two parameters, see the equation (58) of [10], without 
explicit reference to x(q). Here we find it more convenient to introduce, along the line of 
the physical discussion given above, two characteristic "blob" sizes (actually size squared) 
dmm = R 2 — Qk, ^max = R 2 — Qo in terms of which: 

T 1 

-Foo = -— In [27re<i mm ] + — [f(d min ) - /(0) - d min /'(eZ min )] 

+ f^ max) - I*™ Jf^du, (20) 
where d m i n < d max can be found for a given temperature T from the equations 

0£d '»-7H5=r fl2+ ild^ 2 (21) 

Finally, the Parisi order-parameter function, which takes the values between and 1 and 
is the main measure of the ultrametricity in the phase space, has the following shape 

T f"(d) 

x (d) = ""J [yv/( d )p/ 2 ' Vd G [dmm,d max ] . (22) 

where again we found convenient to perform the overall change q — > d = R? — q in 
comparison with notations used in [10]. This function must be now non-increasing, as 
follows from relating its derivative to the probability ir(d) introduced earlier in Eq.(12), 
see the relation Eq.(75). Technically, this property is precisely ensured by negativity of 
the Schwarzian derivative of f(d), see discussions around Eq.(7). 

The above solution is valid for the temperature range < T < T c , where the critical 
(or de Almeida-Thouless) temperature T c is given in terms of the largest blob size d max 
as: 

(dmax) • (23) 

Above this temperature the solution is replica-symmetric, corresponding to a delocalized 
phase for the particle: no particular region dominates the partition function. The corre- 
sponding free energy is given by: 

= -|ln[2^ s ] + -L [f(d s ) - /(0)] - ~ 

where d s satisfies 

d s = R 2 + ^f'(d s ). (24) 

We now consider specifically a correlation functions f(u) of the form (13). In what 
follows we will use the convenient notations z = (2 In i?) _1 and y = z In (u + a 2 ). As noted 
above, our multiscale logarithmic model ensures that $'(y) > and &"(y) > for any 
< y < 1. We also will assume in our analysis below that the function <&"(y) is finite 
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( < <fr"(y) < oo) and differentiable, but later on will relax those conditions. Simple 
differentiation gives: 



/'(«) 



u + a 2 



1 



+ a 2 ) 2 



[$'(y) - , 



(25) 



Our first goal is to find the largest blob size d max from second equation in Eq.(21), and 
then to determine the critical temperature T c . Introduce the scaling variable y max = 
z\n(d max + a 2 ) and using (25), we obtain the following equation determining y max : 



3 (l-J/max)/2 



2 „ ^max/ £ 



a e 



1 + 



z$"(y max ) 



(26) 



Since we are interested in the thermodynamic limit z — ► 0, we can look for a solution 
ymax(z) as a power series of z. One immediately checks that y max {z) = 1 — zln2 + 0{z 2 ). 
This implies that the largest blob size is of the order of the system radius: d max ~ 



R 2 /2 > a 2 for R 



oo. 



Eq. (25) and (23) then yield the critical temperature given in the 



thermodynamic limit, which simply reads: 



T c = yJv{X). 



(27) 



Physically, at T c , the space breaks up blobs of size o(R) and only a finite number of these 
blobs are visited by the particle. However, within each blob, all sites are more or less 
equivalent. Now we can treat along the same lines the first equation in (21) to determine 
the smallest blob size <i m i n for T < T c . It can again be conveniently written in terms of 



the scaling variable y E 



z In (d m i n + a ), such that: 



a 2 e Vrni 



$'(ymin) - Z<^"(y min ). 



(28) 



This equation determines y m i n > for any temperature T < T c and sample size R = 
exp [l/2z\. In the thermodynamic limit z — > 0, it is again natural to look for a solution y m \ n 
as a power series of z, in which we only retain the first two terms: y m - m = + cz + 0(z 2 ). 
Assuming self-consistently that the solution corresponds to > 0, we see that the first 
factor in (28) can be replaced with unity with exponential accuracy. Due to our assumption 
on differentiability of the function $'(y) we expand around y = v*, and after a simple 
calculation find c = 1. This means that d m j n behaves like <i m i n = eR 2u * for R —> oo, where 
satisfies the equation 



(29) 



Since the function &'(y) is monotonously increasing for y > 0, we find that in the limit 
R — > oo (i.e. z — > 0), the equation Eq.(29) must have a unique solution 1 > v*(T) > 
in the range of temperatures \J $'(0) = T m - m < T < T c = \J ^'(l). In this regime, 
^min ^ (imax- Physically, sites within blobs of size d m \ n or smaller are not resolved by the 
particle, which visits all of them more or less equally. 
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Now we can easily find the free energy in the thermodynamic limit z — ► 0. In 
particular, note that Eq.(25) implies that the last term in Eq.(20) can be conveniently 
written as: 

f4iax r _ _ -i 1/2 rill 1 rt/max -. i„ 

I = - (u + a 2 ) 2 /» ' ^~ 2 =— / [<S>\y)-z<S>"(y)] 1/2 cly (30) 

In the temperature range T m j n < T < T c we can substitute here y m - m = (T) + z and 
Vmax = 1 — z In 2, and expand in z up to linear terms. This gives: 

/ = -2\nR C J^y)dy + T c {l + \n2) (31) 

J U-t 

In the same way we evaluate the remaining terms in Eq.(20), and finally find the leading 
and the subleading terms for the equilibrium free energy: = In RT + 6F, where the 
leading term coefficient T is given by: 

-T/T = v*{T) + [HV * ] ~ 2 H0)] + | yffiyjdy, T min < T < T c (32) 

and the correction term is: 

-5F/T = lnv / 2^+l-^ln2-^lna, T min < T < T c . (33) 

For T > T c the solution of (24) in the limit a <C R — > oo is given by: 

and substituting this to (24) we find that the free energy components are given by: 

-T/T=i + m_im, T>Tc , (35) 



-6F/T = \nV2^+Ul + ^ 



T " lin lna. (36) 



J>2 



Finally, the analysis should be reconsidered for T < T m [ n = ^/$>'(0), where d m i n ~ 
a 2 and the particle localizes even on the smallest scales 0(a). Indeed, assuming that 
generically $"(0) < oo it is easy to see from Eq.(29) that « (2T min /<S>" (0))(T -T min ) -> 
as T — > T m ; n . The solution d m i n = eR 2u * is therefore invalid for T < T m j n where z/* stays 
identically zero. Using Eqs. (21) and (25) we find after a straightforward calculation the 
following correct solution to Eq.(28) as z — ► in this range of temperatures: 

dmin = a 2 T _ 0<T<T min , (37) 



13 



showing that T m j n is indeed a delocalisation transition of the REM type, above which the 
minimum blob size becomes much larger than the small scale cut-off a. The leading term 
in the free energy is given by — > limit of (32): 

T = -2^ yf&Wdy, < T < T min , (38) 
Free energy corrections can also be computed, and for all T < T m \ D are given by 

-5F/T = In + |p - § In 2 + (l - 2^p) In a + \ (l + In Jj_) (39) 



1 / T 



2 

In [ 1 



2 

Comparing Eqs. (39) and (33) we see that the correction term is continuous at the 
derealization transition point T = T mm . 

Last but not least, we can determine the thermodynamic limit of the order-parameter 
function x(d) given by Eq.(22), which determines in a precise way how the particle lo- 
calizes on different scales. To leading order in z we find f"'(u) = —2&(y)/u 3 with 
y = zln(n + a 2 ). Introduce again the scaling variable v = 2 1 1 r ^ with d G [eR 2u * , R 2 /2]. 
Denoting the Parisi order-parameter function x(d) expressed in terms of the new variable 
v as X(y) we see that such function assumes the limiting form: 

XH = ji ^ p , V„ €[,.,!]. (40) 

This completes our solution of the problem for the case of continuous function 3>(y). At 
this point it is rather informative to consider the case of a discrete spectrum of exponents 
u, corresponding to K superimposed logarithmic potentials with: 

K 

g 2 {v) P {v)=Y J 9iKv-Vi), Q<v k <v K - 1 <...<v 1 <vq = 1, (41) 
i=i 

with 8{u) standing for the Dirac delta- functions. The corresponding 3>'(y) consists of 
steps: &(y) = Y^f=i 9i@(v ~ u i)- A simple consideration shows that our earlier analysis 
for the values of d max and the critical temperature T c still hold for such a case, so d max = 
R 2 /2, and T c = [^'(l)] 1 / 2 = ^ g\ + g% + . . . g 2 K . The equation (21) used to determine 
^min = R 2ymin — a 2 now takes the following form: 

_ * 1 - a 2 e -y mi n/ 2 i 

~^ 9i l + e (^-J/ min )/,' Z -2l^R- { > 

A little thought shows that the solution should always be in the form y m ; n = v v + c p z for 
small z, where the index p runs successively through the values 1, ...,K when decreasing 



14 



temperature from T c towards T = 0. Introducing a decreasing sequence of characteristic 
temperatures T p = ^Jj2iL p gj, we find the coefficients c p and the index v p for a given 
temperature: 

T 2 - T 2 1 

T p+1 <T<T P (43) 

Thus, the value of y m j n jumps (and thus the size of the smallest frozen blobs <i m j n ) when 
crossing each of the temperatures T p , p = 1, 2, . . . , K with the highest one being T\ = T c . 
Since T p and v p decrease as p increases, it is clear that the order-parameter function X(v) 
for a given temperature T < T c is step-wise constant with jumps at each u p ; the smaller 
v (i.e. the smaller the size of the blobs), the larger X(u), meaning that the condensation 
effect on the scale v is weaker and the Boltzmann weight becomes delocalized for scales 
such that X(v) > 1. Explicitly, in the temperature range T p+ \ < T < T p we find 



X(u) = l + ^_-lJ^-^) + ^__-_J^-^_ 1 ) + ... 
f T T\ T 

... + (---)e(v-v 1 )--e(v-i). (44) 

Invoking the relation Eq.(75) we then see that the probability II(^) for two independent 
particles to end up at an ultrametric separation v is given for the present model by 




... + ^---)5(v-v 1 ) + -S("-l), (45) 

as long as T p+ \ < T < T p , p = 1, 2, . . . K. In particular, in the simplest case of the two- 
scale model K = 2 assuming = < v\ < 1 for the exponents, we obtain the probability 
distribution quoted in Eq.(18). 

The expressions for (y m i n , y max ) suffice to calculate the free energy expression in the 
thermodynamic limit. One easily finds the leading order contribution to be — = InRJ 7 , 
where in the phase with broken replica symmetry we have for p = 1, 2, . . . , K — 1 

2 v x K 

-TIT = v p + - £>i_i -n)Ti + ^ ("i-i - v i) T p+i <T <T p , (46) 

i=l i=p+l 



and finally for p = K 

-TIT = u k + ^ £><-i - Vi) Ti, < T < T K = g K . (47) 



2 K 



T . 



15 



The corresponding replica symmetric expression valid for T > T\ = T c is given by 



1 K 

-TjT = 1 + ^ E(^-i " "0 T ? ■ ( 48 ) 

i=l 

Interestingly, these expressions reproduce exactly, mutatis mutandis the leading order free- 
energy expressions of Derrida's GREM [8] , with a particularly clear interpretation in terms 
of particle localization inside smaller and smaller blobs as the temperature is reduced. 
Corrections to the free energy can also be found, but the corresponding expressions are 
rather cumbersome and are not universal but model-dependent. In the Appendix B we 
provide the explicit free energy expression for any value of R for a different model with 
K— step RSB, which has the same GREM-like thermodynamic limit as the present model. 

We end this section by a comment on the nature of that latter model, which we 
believe deserves separate mentioning. Disordered Hamiltonians usually analysed in spin- 
glass literature give rise to either K = 1 (1RSB), or to K = oo (FRSB) Parisi patterns, see 
for example the results in the framework of the so-called spherical model of spin glasses 
in [33, 34]. As is easy to show, see [35] and also [10], the general class of the models of 
the type (1) includes the spherical model as a special case. It is therefore can be of some 
independent interest to provide an explicit example of a system of this sort which has a 
K— step version of the Parisi hierarchy as an exact solution, for arbitrary K > 1 (see 
[36] for a recent example of a model with a K = 2 RSB solution). In the Appendix B we 
succeed in constructing such an example in the framework of the model (1) of a particle in 
a random potential when N — > oo even for finite values of the sample radius R < oo (this 
case can be looked at as corresponding to a bona fide p-(soft)spin model with a spherical 
constraint). It comes as no surprise that in the appropriately taken thermodynamic limit 
R — > oo this type of models reproduce again the same GREM behaviour as elsewhere in 
the present paper. 

3.3 Multifractality of the Boltzmann-Gibbs measure 

Important information about structure of the Gibbs-Boltzmann equilibrium measure ppir) = 
-zjjfj exp — (3V(r), [3 = y can be extracted from the knowledge of moments 




In the thermodynamic limit of the sample volume Vj, — > oo one expects typically 

m q ~ V~ Tq (50) 

where the set of exponents T q reflects the spatial organization of the Gibbs-Boltzmann 
weights. For example, if the weights are of the same order of magnitude across the sample 
volume, the normalization condition implies locally pp(r) ~ V£ , and a simple power 
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counting predicts the exponents T q = q — 1. In such a situation it is conventional to speak 
about a delocalised measure. The opposite case of a fully localised measure describes the 
situation when essential Gibbs-Boltzmann weights concentrate in the thermodynamic limit 
in a domain with the finite total volume <C Vl — > oo, and are vanishingly small outside 
that domain. This situation is obviously characterized by trivial exponents r q> o = and 
Tq<o = 00 • Finally, in many interesting situations the exponents r q may depend on q 
nonlinearly, and in this case one commonly refers to the multifractality of the measure. 
The Eqs.(49) and (50) imply the following expression for the characteristic exponents r q 
in the general case 

T q = \q\PT(\q\f3)-q(3F(f3) (51) 
relating them to the appropriately normalized free energy of the system: 

m = - lim ^ • (52) 

An alternative way of characterizing multifractality invokes the so-called singularity 
spectrum function /(a). This function characterizes the number dN(a) = da of 

sites in the sample where the local Gibbs-Boltzmann measure scales as pp(v) ~ V£ a in 
the thermodynamic limit. The definition allows to extract the characteristic exponents T q 
as 

T q = - hm i±-!= — . (53) 

H V L ^oo lnV L 

Note, that the restriction of the integration range by the condition f(a) > is necessary 
in order to remove the rare events found in the vanishing number dN(a) — > of sites 
in the thermodynamic limit [13]. This indeed ensures that the extracted exponents T q 
characterize the typical behaviour of the moments. Performing the a— integration by the 
Laplace method one finds that the two ways of characterizing multifractality, by the set 
of exponents T q or by the singularity spectrum /(a), turn out to be simply related by a 
Legendre transform: T q = a*q — /(a*), q = /'(a*). 

The expressions for the free energy calculated in the previous section allow us to extract 
the multifractality exponents r q from Eq. (51-52) and to investigate the corresponding 
singularity spectrum, see the relation Eq.(53). For the systems with a finite number K 
of levels of hierarchy the singularity spectrum f(a) turns out to be piecewise-parabolic, 
generalizing earlier results obtained in the framework of REM-like model with a single 
scale logarithmic correlations [11, 9]. For the case of an infinite hierarchy of scales in 
the landscape , i.e. for a continuous set of exponents is, the behaviour turns out to be 
rather rich and unusual. Below we give a short account of the results for the example of 
continuous function <&'(y) assuming ^'(O) = for simplicity. 

The associated singularity spectrum f(a) calculated via Eq.(53) is positive in an in- 
terval a G (a m i n ,a max ). The positions of the zeros a m i n ,a max of the function f(a) are 
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given by 

ctmin = -PHP) ~ 2/3 £ \f&ty)dy , (54) 



0L r , 



-PF{0)+2P Jny)dy. (55) 



The function f(a) is symmetric with respect to the midpoint of the interval of interest, 
a m = (a m i n + a max ) /2 = —fiT(JS) > 0, where it has the maximum f(a m ) = 1 as expected. 
Close to this maximum, namely, in the subinterval a £ (a_,a+), where the endpoints 
a± = a m ± 2A({3)f- with A((3) = /3 2 ($(1) — $(0)) the singularity spectrum has a simple 
parabolic shape: 

f(a) = 1 - — — (a - a m ) 2 , a_ < a < a+ . (56) 
4A(P) 

In particular, at the boundaries f(a±) = 1 - (3 2 ($(1) - $(0)). Note that in the REM-like 
limit 5>(y) = g 2 y we have a min / max — > and the parabolic behaviour extends to the 

whole interval of positivity of /(a), in full agreement with the results of [11, 9]. 

At the same time for a ^ (a-,a+) the GREM-like model may show a much richer 
multifractal structure manifesting itself, in particular, via a quite unusual shape of the 
singularity spectrum close to the zeros Oi m i n ,<y. max ^ Eq.(54). To illustrate this fact, we 
consider a broad class of functions <J>(y) behaving at small arguments y <C 1 as <E>(y) « 
C 2 y 2s+1 , where s > and the coefficient < C < oo. In particular, in the limiting case 
s -» we are back to the old REM-like model. The behaviour of f(a) in the vicinity of 
the endpoints a m i n or a max is dictated by the behaviour of the multifractality exponents 
at large q, i.e. tu-kx,. As is easy to see from Eq.(51) this asymptotics is in turn extracted 
from the knowledge of the low-temperature behaviour of the free energy: /3^ r (/3)| J g_ >00 . To 
investigate that limit from Eq.(32) we need to know the low-temperature asymptotics of 
the parameter v* related to T via the equation (29). In this way we find that for the 

chosen $(y) the parameter v*(T) behaves as v*(T — > 0) « \ c y , 2a+i ) ' ^^ er simple 
algebra we find with the required accuracy 

-0^)1^*20 £,J&(y)dy-a c (P), (57) 
where we introduced the short-hand notation 

Now the relation Eq.(51) immediately yields the required asymptotic behaviour of the 
multifractal exponents: 

T f ga min - a c (/3)q~°, for g -> oo 
1 qot-max ~ u c (0)\q\~s , for g -> -oo 
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where a min / max are precisely the zeroes of f(a) given by Eq. (54). We thus see that the 
singular behaviour revealed by (59) is very different from the situation typical for other 
types of disordered systems[13, 12] where one always observes a precise linear behaviour 
T q = ciot m i n / max starting from some value of \q\ (see formula (2.42) in [13] and discussions 
around it). In particular, such anomalous asymptotics is translated by the Legendre 
transform to the anomalous singularity spectrum behaviour close to the left and right 
zero: 

/(<*) ~ s s/(s+i) a c{(3) s,(s+1) \a - a min/max \^r , (60) 
We see that for any s > the derivative of the singularity spectrum diverges when ap- 

s 

proaching zeroes as / (a) ~ \a — ct min / max \ "+ 1 — > oo. Note that this is again very 
different from the standard behaviour observed in other disordered systems where always 
f'(a) < oo at zeros of f(a) [13]. 

4 Discussion and Conclusion 

In this last section, we want to discuss several aspects of our model which are certainly 
worth investigating further, among which: (i) the case of finite dimensions; (ii) the dy- 
namics in such a multiscale landscape and (hi) the relation with multifractal random 
walks. 

4.1 Multiscale logarithmic landscape in finite dimensions 

As detailed at the end of Sect. 2.2, the exact, N — > oo results for the single scale loga- 
rithmic potential match perfectly with the results obtained by Carpentier and Le Doussal 
using RG, numerical and heuristic methods [9] in finite dimensions. There is no reason to 
doubt that these results are in fact exact in all dimensions N > 1. Although our model is 
substantially more involved, we see that essentially the same physical mechanisms are at 
play in both models, in particular in the case of a finite number K of hierarchies. There- 
fore, it is very tempting to conjecture that the GREM behaviour revealed by above in the 
infinite-dimensional setting should also hold in all spatial dimensions, down to N = 1. It 
would be very interesting to see if the corresponding RG and travelling wave formalism 
can be generalized to support this conclusion. In the case of finite K, this looks indeed 
quite feasible. 

If this conjecture is true, we would then have indeed explicitly constructed a Parisi 
landscape in finite dimensions fully in terms of stationary Gaussian processes. How do 
we reconcile this with the ultrametric properties of the Parisi construction? Consider the 
following distance Dr defined for any two points r, r' inside a sphere of the radius R in 
the Euclidean space of any dimension N: 

g,(ry)- ' n P r ^ +a 1 , 0<|r|,|r<|<iJ (61) 
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Parameterizing |r| = R a ( r \ < a < 1, we see that in fact 

lim D R (r, r) = max{a(r), a(r')} . (62) 

R— >oo 

One can easily check that the latter function used as a distance converts the Euclidean 
sphere into a so-called ultrametric space: every triangle will have at least two sides equal. 
We thus conclude that our choice of the model corresponds to insisting that the covariance 
of the random potential values in two points in space should depend only on the ultrametric 
distance inside our growing sphere. The original construct of the GREM by Derrida in 
fact proceeds in a similar way, and therefore the coincidence between our model of the 
GREM could have been, with hindsight, anticipated. Indeed, Derrida started with 2 M 
random variables, attached to the vertices of a hypercube of 2 M points supplied with a 
tree structure and associated ultrametric distance [31]. The Gaussian random energies 
used for constructing the Boltzmann weights were built from those ingredients in such a 
way that their correlation function depended only on that ultrametric distance. Although 
the methods used in this paper have very little in common with techniques used by Derrida, 
or other authors analysing GREM, the similarities between the two problems is apparent. 
This line of reasoning suggests that the convergence to the GREM limit could be indeed 
expected by invoking universality arguments. A remarkable recent progress [31] based on 
the idea of Ruelle probability cascades [37] has allowed one to analyze Derrida's original 
GREM in full mathematical rigour. It would be extremely interesting to see whether the 
Euclidean version of the GREM analysed in the present paper is amenable to a similar 
kind of rigorous analysis, without any reference to the replica trick, powerful heuristically 
but still ill-defined mathematically. 

Another, more speculative aspect of the problem is also worth mentioning here. In a 
recent work, Moore [38], argued that low-temperature phases with 1RSB fail to survive 
in finite spatial dimensions. This, in Moore's argument, is intimately related to a generic 
absence of marginally stable modes in a 1RSB fluctuation spectrum. To this end, the sta- 
bility of the 1RSB low-temperature phase in the present model (1) was investigated in [10], 
and found to be controlled by two eigenmodes, denoted in [10] as Aq and A^ (see equations 
(B.29) and (B.30) of the Appendix B of [10]). Generically those two were demonstrated 
to be positive, but for the logarithmic potential (8) both Aq and A^ identically vanish ev- 
erywhere in the low-temperature phase. That property makes the associated 1RSB phase 
marginally stable 7 . Thus, the glass behaviour in the case of logarithmic correlations could 
survive for finite N due to marginality of the fluctuation spectrum. This picture would be 
indeed in agreement with the above-discussed RG results of [9]. We believe that similar 
marginality of fluctuations spectra should exist also in the general K— step case of our 
model, and is related to the fact of vanishing Schwarzian derivative, exactly or asymptot- 
ically in the thermodynamic limit. The corresponding calculation looks cumbersome, but 

7 This fact, though not explicitly mentioned in [10], immediately follows from definitions (B.29) 
and (B.30) after substituting for qi — qo = Q and qd — qi = y the expressions (74) and (79) of that 
paper. 
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is certainly feasible, see similar work in [39], and is yet to be done. In general, any results 
in these directions are highly desirable, with an ultimate goal of performing perturbative 
expansions around N = cc limit which remains to be one of outstanding challenging tasks. 

4.2 Diffusion in a multiscale logarithmic landscape 

The rich behaviour found in the thermodynamics of a single particle in a random potential 
also has interesting dynamical counterparts. In the infinite dimension limit, the problem 
of a Langevin particle in a random potential has been solved in details by Franz and 
Mezard [40] and Cugliandolo and Le Doussal [35], for both the short-range and the long 
range cases. These results reveal long-time relaxation, aging, and other effects typical of 
glassy dynamics. The marginal (monoscale) logarithmic case was however not specifically 
studied in these papers. For finite N, this marginal case was studied in [42, 43] using RG 
methods. The main quantity of interest is the dynamical exponent z, defined as: 

A 2 (i) = <(r(t + t ) - r(t )) 2 ) ~ t ^oc t 2 ' z . (63) 

The result, conjectured to hold at all orders in perturbation theory and for any N, is that 
z varies continuously with the strength of the disorder g and temperature T: 

z = 2 + 2 (|) 2 . (64) 

When g = 0, one recovers the diffusion exponent z = 2, as it should. It was noted 
by Castillo and Le Doussal [41] that this result cannot hold down to T = in N = 1 
dimension, because it violates exact bounds. It was argued instead that below the static 
transition T c = g found in Sect. 2.2, the exponent in fact is modified and reads z = 4T C /T. 
Therefore the static transition is also a dynamical transition in N = 1. 

The situation in higher dimensions is unclear, in particular it is not immediately clear 
that the dynamical transition where Eq. (64) ceases to hold still coincides with the static 
transition. This is another motivation for studying in details the limit N — > oo. However, 
it seems reasonable to conjecture that for the multiscale landscape model at high enough 
temperature, the diffusion exponent becomes scale dependent. For a finite number K of 
hierarchies, and for T > T\ = T c , we expect to find: 

z(A) = 2 + 2^) 2 , (65) 

where T p is the critical temperature associated to u p such that R Up < A < R^p- 1 . At 
temperatures lower than T c , it would appear natural to conjecture that all levels such that 
T < T p follow the Castillo-Le Doussal scenario, whereas levels such that T p > T are still 
ruled by the above exponent. However, the situation might be more complicated because 
the possibility of aging at low temperatures seems to have been overlooked in [41]. In high 
enough dimensions, the analogy with the GREM suggests that the multiscale logarithmic 
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model of the present paper might be a real space realisation of the multi-level trap model 
introduced and studied in [44]. We would then witness a very rich dynamical behaviour, 
where all levels such that T < T p age (concerning large length scales) , whereas small length 
scales, such that T > T p , are still stationary [44, 45]. It would be extremely interesting to 
study these aspects in more details. 



4.3 A Generalized Multifractal Random Walk 

The monoscale logarithmic landscape model in N = 1 has in fact deep connections with 
the muZMractal Random Walk (MRW) construction of Bacry, Muzy & Delour [14, 46]. If 
one treats the coordinate r as a time variable t, and call X t the position of a random walk 
at time t, the BMD model is defined by the following evolution equation: 

dX t = mdt + atdWt at = ao exp f3Vt, (66) 

where dW is the usual Brownian process, m a drift and Vt is the logarithmically correlated 
Gaussian process considered above (possibly shifted such that (e 2l3Vt )v = !)• It was shown 
by BMD that the resulting process is multifractal in the time regime a <C t <C R. For 
m = 0, all odd moments vanish and even moments are given by [14]: 

{(X t -X t+T r) = M n r^ with Cn = ^-(/3<7) 2 ^ 1 ^, (67) 

as long as n < ([3g)~ 2 beyond which all moments diverge: the distribution of increments 
A = (X t — X t + T ) has a power-law tail A^ 1 ^^ with an exponent given by: 

,= (|) 2 . (68) 

The connection with the thermodynamical problem addressed above is clear: the MRW is 
a standard random walk subordinated to a stochastic time St, such that: 



rt+r 

■H+r-st = J due 2pv \ (69) 



which is a restricted partition function at temperature T/2. 

The generalisation to a multiscale logarithmic process is very natural. Taking Vt = 
Y^v=i Vv,t as i n Eq. (14) above, one constructs a Generalized MRW with, in the limit 
R — > oo, K different epochs such that v p < lnr/lni? < v p -i such that the moments of the 
walk are given by the above multifractal scaling, Eq. (67), but with a different multifractal 
spectrum: 

C« = f-(f) 2 ^. (70) 

This formula can be given a meaning in the continuous limit as well. In this case the 
exponents (n^ acquire a continuous dependence on lnr: they are still given by Eq.(70) 
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but with the ratio T/T p replaced with the Parisi function X(v), v = lnr/lni?. It would 
be interesting to understand in more details the extreme value statistics of the GMRW, 
following the method introduced in [9, 47]. 

4.4 Final remarks 

The model we have introduced appears to contain a very rich phenomenology, and deserves 
in our opinion further investigations. Many points would require a more rigorous mathe- 
matical investigation, concerning in particular the finite dimensional version of the model 
and its dynamical properties, in particular in the continuous hierarchy case K — > oo. We 
believe that our construction sheds an important light on the understanding of replica 
symmetry breaking: clearly the fact that Parisi's ultrametric scheme is relevant for a 
translationally invariant, Gaussian potential in finite (even one) Euclidean dimension is 
very satisfying, if only for pedagogical reasons. A better intuition on the Parisi solution 
for the SK model for a large but finite number of spins M, and the recently reported 
scaling K ~ M 1 / 6 [48], might also be within reach. 
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Appendix A: Calculation of the overlap function, 
Eq.(12) 

In the framework of the replica trick we represent the normalization factor in the product 
of the Boltzmann-Gibbs weights as l/Z 2 = lim n ^oZ n ~ 2 . This trick allows the disorder 
average to be performed explicitly, and (12) takes the form 

/I \ n 

^ / (2=1 
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a=l 

At the next step we exploit the 0(N) rotational invariance of the integrand and, 
assuming N > n, introduce the scalar products r a rfe = q a b as new integration variables, 
see [10] for a general description of the method. After further rescaling q a b — > Nq a b the 
integral in the right-hand side of the equation takes the form 

C N , n N N ^ J d (detQ)-^/ 2 e~ N ^Q) 5 (v _ * {qu + , 22 ) + ^ d Q (71) 



where 



$ n (Q) = -^ln(detQ)-/3^/ 

' a<b 



-^{qaa + Qbb) - Qab 



(72) 



and the integration domain Dq over the matrix Q with entries q a b is already N— independent: 
Dq = {Q > 0, q aa < R 2 , a = 1, . . . n}. The proportionality constant Cn,u is also known. 

The form of the integrand in Eq.(71) is precisely the one required for the possibility of 
evaluating the replicated partition function in the limit iV — > oo by the Laplace method. 
Taking into account the permutational symmetry of the integrand with respect to the 
replica indices, we find after a straightforward calculation 

tt(P) = lim — — Yl 6 ( V ~ 2^ qaa + %b) + q A iQestationary point ■ ( 73 ) 

To perform the replica limit explicitly we rely upon the assumption of validity of the 
hierarchical ansatz for the stationary-point solution suggested by Parisi, see a description 
in the context of the present model in Appendix A of [10]. Denoting mi the size of the 
blocks in the Parisi scheme, we find in the standard way (cf. Eq.(46) of [10]): 

tt(V) = lim V(m /+1 - mi)5(T> -q d + qi )= / 5{V -q d + q)x'(q) dq = x'(q d - V) , (74) 

where x(q) is a non-trivial non-decreasing Parisi function characterizing the low-temperature 
phase. In the present model it turns out that the parameter q d is related to the effective ra- 
dius of the sample as q d = R 2 - It is actually more convenient to use the variable d = R 2 — q 
rather than q itself, as we do in the main body of the paper. Accordingly, we pass from 
x(q) to the non-increasing function x(d) to which tt(V) is especially simply related: 

n{V) = ~x(p), (75) 

which is the counter-part of the standard interpretation of x(q) in terms of the overlap 
probability in spin-glasses ([49]). 
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The above construction is valid for any choice of random Gaussian potential with any 
covariance function of the form Eq.(l). For the case of multiscale logarithmic potential 
studied in the present paper the argument d of the function x(d) generically covers the 
interval d £ [d m i n , d max ], with d m i n ~ R 2u * and d max ~ R 2 , where the value determined 
for a given temperature T by Eq.(29). Accordingly, it is natural to pass from the Parisi 
function x(d) to its counterpart X(u), v E by replacing d — > R 2v . As is clear, v 

is precisely the " ultrametric" separation between the two positions related to the scaled 
squared Euclidean distance T> as v = 2 \nR In^l-R— >oo) cf. Eq.(62). It is then immediate to 
check that the probability Il(^) for two independent particles to end up at an ultrametric 
separation v is related to the function X(y) in precisely the same way as 7r(£>) to x(d): 

U(u) = -2\nR^L- )X {V)\ v=R ^ = ~X{y). (76) 

Appendix B: A model for K-step replica symmetry 
breaking 

The idea of the construction comes from an observation that the most general function 
f(u) with locally vanishing Schwarzian derivative S(u) = 0, see Eq.(7), is given by f(u) = 
f — Au — g 2 \n(u + a 2 ). Apparently, linear in u term violates the requirement f'(u — ► 
oo) — > and for this reason was discarded by us earlier. Now we can try, by choosing 
constants / and A judiciously, to construct a covariance function f(u) globally from the 
above local patches in such a way, that the composite solution will have the vanishing 
S(u) everywhere, except in a finite number of points. As the derivative x'(d) of the Parisi 
order-parameter function from Eq.(22) is proportional to the Schwarzian derivative S(d), 
this construction should provide us with a piecewise-constant order-parameter function 
required for the Parisi pattern with K levels of hierarchy. For the sake of simplicity we 
discuss below in detail the simplest non-trivial case K = 2. The generalization to arbitrary 
K > 2 will be apparent. 

We find it more convenient to work with the structure function 4>(u) = /(0) — f(u) 
rather than with f(u) itself. The structure function obviously vanishes at the origin, pro- 
viding an additional condition 0(0) = used to specify all the constants in our construction 
uniquely. 

Consider the model of a particle in a Gaussian random potential with the structure 
function taken in the form 

J /i- A lU - t 2 In (u + aj), u>u* 
n> ~\f2-A 2 u-t 2 ln(u + a 2 ), < u < u* ' 1 > 

where the positive parameters t±, a±, a 2 are considered to be given, and chosen satisfying 
the inequalities: 

a 2 a 2 

ti > t 2 and -i > -2. (78) 
h t 2 



25 



In contrast, the value of the "breakpoint" «*, as well of the constants f\, / 2 , Ai, A2 are 
unknown and should be specified in terms of 01,2 and ii j2 . Note that for <p{u) to have the 
meaning of a structure function requires A\^ to be non-negative, see eq. (3). 

As we are to use such a function for building the Parisi order-parameter x(d) according 
to Eq.(22), and such x(d) can take only finite values between and 1, one has to ensure 
that 4>(u) is continuous together with at least two derivatives everywhere, including the 
breakpoint -u*. The requirement of continuity of 4>"(u) at u* immediately fixes the value 
of the breakpoint: 

o2 _ o| 

u * = TTT' (79) 

*2 tl 

The consistency of the procedure requires -u* > which is precisely the case due to the 
imposed conditions Eq.(78). 

In the same way the continuity of 4>'(u) at given by (79) allows us to relate A\ to 
A2 via 

A 2 = A 1 + %^ (80) 
af — a 2 

At the next step the continuity of (j)(u) at u* relates f\ to / 2 as 

h = h + i h ' (tl - *2) + tl lntx - t| lnt 2 + (t? - «2) In . (81) 

The constant / 2 can be found from the condition (f)(0) = 0, which fixes / 2 = t 2 lna 2 . 
Finally one should have in mind that although for a finite sample x < R 2 , to have a well- 
defined model for any x < 00 we should impose the condition <p'(u — > 00) = as elsewhere 
in the paper. This necessarily implies A\ = 0, and fixes all the constants of the model 
uniquely. Note that A 2 is indeed positive in view of Eq.(80) and the inequality a\ > a 2 . 

Having specified the structure function <fi(u), hence the covariance f(u), we now can 
repeat the procedure of finding <i rmn , d max , and hence the free energy of the model. As- 

R 2 —a 2 

suming d ma , x > u* we can solve the second of equations Eq.(21) and find d max = — 2 - Ji - 
Such solution is indeed the valid one for sufficiently large sample sizes satisfying R 2 > 
= [ a i(^i + £2) — 2a 2 ti]/(ti — i 2 ). Assuming the condition is satisfied, the transition 
to the phase with broken replica symmetry occurs at the temperature found from Eq.(23) 
which is given by 

?2 „2 



R z -gj 
R 2 + a 2 

Then the solution to the first of equations Eq.(21) reads: 



rc = tl-S5-TT3- ( 82 ) 




<in= ' 2 %V mm T' , (83) 
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and further requiring consistency of this solution with the expression Eq.(79) finally yields: 

"5 f / 1 1 - f ) • T »» £ T < T = 

where the "breakpoint temperature" T m j n is given by 

T m i n = t\ t2 o • (^5) 

It is easy to see that indeed £2 < Tmin < T c < t\ for R > i? m , so such a solution is 
consistent. Note also that d mm (T) < u* for T < T m i n , and d m i n (T c ) = d max as expected. 

A simple calculation shows that the Parisi order-parameter function x(d) in the tem- 
perature range < T < T rmn indeed consists of two perfect steps: 

/ j\ / J for d min < d < U*, 

»(rf) = i T n „. ^ ,^ , • (86) 



tl 



for u* < (f < d max 



whereas in the range T imn < T < T c the value d m i n exceeds the breakpoint u*, and hence 
only a single step survives: 

T 

x(q) = — for d min < d < d max . (87) 

Now we can easily calculate the free energy for T < T c from Eq.(20). In particular, as 
we have for < T < T min 



yfrvi-\ % ,„<„<„ < 88 > 



we find that in this range of temperatures 

/ = - /*~ [/»] 1/2 d« = -fc In f -^U - h In f *=!S±£ 
At the same time, for T m i n < T < T c we have y/ f"(u) = so that 

\ Omin + af J 

The free energy Eq.(20) is then given by 

T 2T 2 2 ti T 



2 2 



(89) 



2 v t ; i-f 



T < T < T 

± mm - 1 ^ J c 
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and 




h 1, T t 2 
— k H — In 

2T 2 2 t 2 T 



- In (u* + a|) 



T 




- - In (u* + a\) 
< T < T min (90) 



One can check that at the breakpoint temperature T = T rmn the free energy is continuous. 

Finally, for T > T c the replica-symmetric expression for d s obtained from Eq.(24) is 
given by 

R 2 - a 2 + ^(i? 2 + a 2 ) 2 + AR 2 a 2 t 2 /T 2 



(L 



2(l + i 2 /T 2 ) 



(91) 



and the free energy is obtained by substituting this to Eq.(24), and remembering d s > «*■ 
One can check that cI s \t^t c = (R 2 — «i)/2 = d mSLX , and the free energy is again continuous 
at the transition temperature T c . 

All the above expressions were derived for finite sample sizes, provided R ^> Rm- 
Let us now investigate for this model the thermodynamic limit R — > oo, taken at a fixed 
temperature T. As elsewhere in the paper, we assume the scaling a\ = R U1 ,a 2 = R U2 , with 
< v 2 < v\ < 1, which implies ai 3> 02 in the thermodynamic limit. We immediately find 
that the critical temperature tends to T c = t\, and d max — ► R 2 /2, d s — ► R 2 /(l + T 2 /T 2 ). 
Similarly, the breakpoint temperature T m i n — > t2 5 and the breakpoint itself is given by the 
limiting value u* = i? 2 " 1 /(T c /T mm — 1). We also need to know that = 2f 2 T^ in lni? and 
/1 = 2 [z/i(T 2 - T^ in ) + ^ 2 T r 2 in ] lni?. The leading contributions to the free energy found 
from Eqs.(89) and (90) are given by 



TlnR 



~^( Ul - U2 ) + ^(1 _ < t < r min 

T 2 



"2 + — 

V X + ^(1 - I/x) + (U! - U 2 )^, 



1 < T < T 

mm ^ ^ J-c 



1 + (1-^1)^ + ^1 -"2)'%", T>T, 



rp2 



(92) 



These expressions are identical to K = 2 case of GREM free energy (46) and (48), upon 
identification T m - m = T 2 ,T C = T\ , showing that the two models are indeed in the same 
universality class in the thermodynamic limit. The corrections are however already model- 
dependent. 

We further note that the parameter d m i n given by Eq. (83) develops in the limit R — > 00 
characteristic divergencies at the breakdown temperature T m i n as well as at T c (cf. (43)): 



(93) 



r_ R 2v 2/ f 1 _ y T_\ 0<T<T mil 

am ' \ J mm / 

^i^V(l-f), T min <T<T c 



The corrections to the free energy however do not have any divergencies as the logarithmic 
terms come with vanishing prefactors. 
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The K = 2 construction discussed above has obvious generalization to any K > 2. We 
choose the structure function in the form 

-<p{u) = f p - A p u-t 2 p \n{u + a 2 p ), < u < ul P ~ 1} , p = 1, 2, . . . , K (94) 

where the positive parameters t±, . . . , tx', 0,1, ■ ■ ■ , clk are given, and chosen satisfying the 
inequalities: 

h>t 2 >...>t k and ^>^>...>^. (95) 

Those inequalities will ensure that a strictly decreasing sequence of the breakpoints = 
R 2 > > > . . . > ui K ^ > can be found from the conditions of continuity of 
(f)"(u). Requiring also continuity of (j)'{u) and <j)(u) relates A p to A p _\ and f p to f p -i for 
any p. Finally we set </>(0) = which gives fx = tj^lnaj^, and put A± = to satisfy 
the behaviour at u — > 00. This choice fixes all the parameters of the model uniquely, and 
ensures A p > 0. Considering the effective radius R exceeding some minimal value R m ~ a 2 , 
the model will have the sequence of "break temperatures" T c = T\ > T2 > . . . > Tk > 0, 
with the order-parameter function consisting of / decreasing steps for < T < 7), with 
the values at those steps given by x{d) = T/ti, for u$ < d < ui ^ i = 1, 2, . . . , I, where 
we make a convention = d max , d[ 1 ^ = d m i n . Taking the thermodynamic limit R — ► 00 
with the the powerlaw scaling of the cutoffs as = R 2ui reproduces exactly the GREM 
equations (46) and (48). 
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